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DAMPING OF A TORSIONAL OSCILLATOR IN LIQUID 
HELIUM 4 AND 3 FROM 0. 4° TO 2. 5° K 

by Rayjor W. H. Webeler and David C. Hammer 
Lewis Research Center 

SUMMARY 

The logarithmic decrements of a piezoelectric, cylindrical, quartz crystal driven in 
a torsional mode of vibration were determined as a function of temperature with the 
crystal immersed in liquid helium 3 (He^) and in helium 4 (He^). These measurements 
were then used to calculate the applicable "viscosity -density" product of each of these 
two liquids as a function of temperature. 

The results of this investigation in the interval 2. 14° to 2.20° K represent the most 
extensive and precise viscosity-density values for He that are available to date. Inspec- 
tion of this viscosity -normal density product (T? n P n ) data near the lambda point (T^) in- 
dicates that a discontinuity in r? may exist. 

Results from this investigation between 2. 0° and 1. 35° K in liquid He* are in reason- 
ably good agreement with three other investigations. Between 1. 1° and 0. 95° K, the data 
presented here substantiate to a large extent only the rotating viscometer results of 
Woods and Hollis Hallett. At temperatures below 0. 6° K where there are no comparable 
investigations, this crystal decrement has a different significance than at higher temper- 
atures. In this temperature region this investigation found the phonon density temper- 
ature dependence to be T^' confirming the temperature dependence of the Landau 

expression for the phonon density. 

No evidence was found that the torsional crystal technique provided a means of 
separating the phonon and roton viscous contributions, but there was evidence of suppres- 
sion of the phonon viscous contribution. 

Viscosity values for liquid He^ were determined between 0. 4° and 2. 5° K and are 
compared with the results of other investigators. Below 0. 6° K, these data are in ex- 
cellent agreement with those of Betts, et al. , but they disagree at higher temperatures. 


INTRODUCTION 


In terms of a simple working model, liquid helium 4 (He^) below the X point con- 
sists of a normal component and a superfluid component. The total density of the liquid 
is the sum of the superfluid density and the normal density. The superfluid component 
does not take part in the momentum transfer or viscous properties of the liquid so that 
only the normal component has a nonzero viscosity coefficient. The Landau theory 
applied to this model would consider the normal component to consist of two types of 
thermal excitations called phonons and rotons which transfer momentum as a result of 
various scattering processes. The normal component may be described as consisting of 
"particles” with rotons represented by relatively heavy particles and phonons by rela- 
tively light particles. The density of both rotons and phonons is highly temperature 
dependent. The roton density is higher above 1° K than the phonon density. Below 
0. 6° K, the normal component consists mostly of phonons. Thus well below 0. 6° K, 
phonon-phonon scattering predominates. 

Landau and Khalatnikov (see ref. 1, pp. 106-110) have considered the normal com- 
ponent as a gas of phonons and rotons with very weak interactions below 1. 6° K so that an 
ideal gas model may be used. The phonons and rotons are treated as particles. Landau 
and Khalatnikov divide the viscosity of the normal component into a phonon viscosity 

and a roton viscosity Assuming that phonon-roton collisions ensure the same 

drift velocity for both phonons and rotons, 77 n = ? 7 R + T 7 ph where r/ n is the total viscosity 
of the normal component. The roton viscosity is due primarily to roton-roton collisions. 
The phonon viscosity is due to phonon-roton and phonon-phonon collisions. In addition 
they have suggested that the roton viscosity is plausibly temperature independent. The 
phonon viscosity on the other hand is highly temperature dependent. It may be expressed 
in the form = ap^c Z p ^ where a is a constant of proportionality, p p ^ is the phonon 
density which decreases rapidly with decreasing temperature, c is the velocity of ordin- 
ary sound in the liquid, and Z p ^ is the phonon mean free path which increases rapidly 
with decreasing temperature. Above 0. 9° K, Z p ^ is determined primarily from phonon- 
roton collisions. At lower temperatures where the roton density decreases rapidly, 
phonon-phonon collisions determine Z^. With decreasing temperature, Z p ^ increases 
more rapidly than the phonon density decreases, so that the phonon viscosity increases. 

4 

If the roton viscosity is constant with temperature, the total viscosity for liquid He will 
be nearly temperature independent just below 1. 8° K and then increase as the temper- 
ature decreases, which agrees with experimental results. 

The present work was motivated partly by the discrepancy in viscosity data below 
1.0° K between a mechanical rotation technique and one which measured the attenuation 
of second sound from viscous effects at the wall of a resonant cavity. Also Atkins 
(ref. 1, p. 110) suggested the torsional crystal technique may partially suppress the 
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phonon viscous decrement but measure the full roton viscous decrement. This would 
provide a method of separating the roton and phonon viscous decrement. This sugges- 
tion further stimulated interest in the temperature region just below 1° K, where this 
effect would be most evident. 

4 

The A. point transition in liquid He can be carefully examined with the torsional 
crystal technique as demonstrated by Welber and Quimby (ref. 2). We believed that even 
better resolution was possible in this region giving very precise viscosity -density mea- 
surements in the temperature interval on either side of the transition. 

Liquid He^ is a system of bosons which may be considered to undergo an imperfect 

Q 

bose condensation at its A. point while liquid helium 3 (He ) is a system of fermions 

3 

which from first principles is not expected to undergo such a condensation. Liquid He 
is an isotropic collection of like particles having no charge and spin 1/2 with strong 
interatomic forces. The Landau theory for a strongly interacting system of fermions at 
temperatures below 0. 1° K predicts a viscosity coefficient which varies as T - ^. Mea- 

n 

surements of the viscosity coefficient of liquid He J below 1. 0° K are available from two 
previous investigations, but the values disagree by 20 percent. The present investiga- 
tion was undertaken to resolve this discrepancy. 


SYMBOLS 


C’ additional capacitance 

C capacitance in parallel with crystal 

px 

C_ capacitance 

C x capacitance 

Cgg applicable rigidity modulus 

c velocity of ordinary sound in liquid 

E root mean square voltage across crystal at resonance, IR, V 

F resonant frequency of ideal torsional rod 

F g electrical resonant frequency of series branch 

Fj^ mechanical resonant frequency 

I root mean square current, A 

K q "constant" of crystal system 

2L total length of rod 

L x inductance 
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ph 


N 

P 

P 

R e 

R s 
R. . 


R 


X, vac 


S 

T. 


M 

V m 
v(x, t) 

V S 

W d 

W v 

z,z a , 

Z b’ Z n’ Z s 


a 

6F 

A 

A. 

i 

V 


2R 


% 

% 


phonon mean free path 
mass of crystal, M 

3 

number of particles per cm 
classical momentum per particle 
pressure 

radial coordinate of end faces 
variable bridge resistance 

measured crystal resistance at resonance, ohm 

crystal resistance in vacuum corresponding to temperature at which R is 

X 

measured 

radius of crystal cylinder 
total surface area of crystal 
temperature at X point 

maximum velocity of point on rim of cylinder 
maximum velocity of point on element of area 

velocity of particle in fluid at distance x from plane surface at time t 

velocity of shear in quartz 

energy dissipated per cycle, J/cycle 

energy of vibration (maximum) 

impedance 

constant of proportionality 

width of resonance curve at half power point 

logarithmic decrement 

logarithmic decrement due to fluid surrounding crystal 

viscosity coefficient of liquid (shearing stress/velocity gradient normal 
to flow) 

total viscosity of normal component, 7/ n = *7 ^ + 77^ 
phonon viscosity 
roton viscosity 
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9 polar angle 

9 instantaneous angular velocity 

Qjyj- maximum angular velocity corresponding to point on end face or rim of cylinder 
1/A C. B. L. thickness (depth) 
p density of liquid (fluid) 

P n density of normal component of liquid helium 4 
Pq density of quartz crystal 
Ppk phonon density 
p n roton density 


ANALYTICAL CONSIDERATION 

A cylindrical alpha quartz crystal may be driven in a fundamental torsional mode of 
vibration by the proper application of an alternating electrical field. Such a crystal, 
ground in the shape of a right circular cylinder, must have one of the three "electric 
axes" of its lattice alined parallel to the geometric cylindrical axis in order to produce 
a fundamental torsional mode of vibration. This torsional mode will be relatively pure 
provided additional care is taken to aline the alternating electric field with respect to the 
"optic", "x", and "y” axes of the quartz crystal. Such alinement of the alternating 
field with respect to the crystal lattice is discussed by Mason (ref. 3, p. 91) in his 
treatise on piezoelectricity and need not be repeated here. 

The fundamental torsional mode of vibration consists of a strain node at the midplane 
perpendicular to the longitudinal axis of the quartz cylinder with a stress node at either 
end. The crystal is suspended freely at the central strain node. 

When this cylindrical quartz crystal is oscillating at resonance in a fluid such as 
liquid helium, energy is dissipated because of 

(1) Internal friction of the quartz crystal, 

(2) Mounting or suspension losses, and 

(3) Viscous forces of the liquid 

Additional sources of energy dissipation arise from sound waves because of imperfect 
alinement of the crystal and from deviations in the shape of the crystal from a perfect 
right circular cylinder during torsion (ref. 4). Specifically when half the circumference 
of the crystal or half the inside circumference of the electrode boundary wall is not less 
than one wavelength of sound in the sample fluid, an unwanted resonant condition will 
probably exist at certain sample fluid temperatures. This will result from the formation 
of standing sound waves around the inside of the boundary wall, with one or more integral 
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number of wavelengths around one-half the ’'crystal" circumference. This resonance 
may be eliminated by selecting a crystal with a satisfactory small ratio of diameter to 
length and placing the electrode boundary wall near the crystal surface. 

As a consequence of the previously mentioned causes of dissipation, electrical energy 
must be supplied externally in order to maintain a constant amplitude at resonance. One 
expects the crystal amplitude and the power dissipated to be a maximum at resonance for 
a given peak alternating voltage of the driving electric field. The "viscosity times 
density" of the surrounding fluid is a function of the logarithmic decrement A for the 
crystal system. 

The logarithmic decrement of a piezoelectric crystal or of any mechanical vibratory 
system is defined by the quotient of energy dissipated per cycle divided by twice the 
energy of vibration W . In this case the energy of vibration is simply the maximum 
kinetic energy attained by the crystal during a cycle. 

Mechanically this crystal is simply a torsional rod clamped at the center with both 
ends free such that these ends may twist in opposite directions. If the total length of 
such a rod is taken as 2L and this rod is oriented such that the applicable rigidity 
modulus is Cgg, then the resonant frequency of the torsional rod F is given by (ref. 3, 
p. 67) 



where Pq is the density of the quartz crystal. 

The mechanical resonant frequency F R for this crystal rod will be nearly equal to 
the resonant frequency of the ideal crystal described. The equivalent circuit for the 
crystal in the system (fig. 1) consists of a series resonant branch having an inductance 
L , capacitance C , and resistance R ; these being in parallel with an additional 
capacitance C'. (A detailed analysis of crystal equivalent circuits is given by both 
Mason and Cady, refs. 3 (p. 67) and 5, respectively). The electrical resonant frequency 
for the series branch of this circuit is well-known and given by F g = (1/2it)(1/VLC). 

This is also equal to the resonant frequency of the entire circuit provided it is defined 
as that frequency at which the external current flow is a maximum. There is also a 
frequency for the equivalent circuit at which the external current flow is a minimum. 

This is called the antiresonant frequency. In a high "Q" system such as this one, this 
antiresonant frequency will be very near to but slightly higher than the resonant fre- 
quency F r . 

At the resonant frequency of the series branch, the crystal exhibits only electrical 
resistance and no reactance. Under these conditions there is a relation between viscosity 
times density of the surrounding fluid and the crystal resistance in ohms (see appendix). 
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(b) AC capacitance bridge for viscosity times density measurements. 
C" - C px + C\ X x represents crystal reactance (0 at resonance), 
and Z a , Zjg, Z s , Z u are arm impedances. 

Figure 1. - Equivalent circuitry. 


7 


II I I I 


I I 


I 


I llll I I « I 


III ■■lllll■■■l 



This relation is 

VP = 

K o S V 

which is a special form of the more general equation 

a 2 = jr&ip (A10) 

t P 

F R M 

where is the logarithmic decrement of the crystal due to the surrounding liquid, S 

is the total surface area of the crystal, 7j is the viscosity coefficient of the liquid with 

density p, F-^ is the resonant frequency of the crystal system, and M is the crystal 

mass. R x is the measured crystal resistance at resonance in ohms with the crystal 

immersed in the liquid of interest, and R is the crystal resistance in vacuum cor- 

x, vac 

responding to the temperature at which R„ is measured. R v corrects for dissipa- 

x x 5 vac 

tion in the crystal due to internal friction and suspension losses. K q is a "constant" 
of the crystal, and its determination will be considered later. 

In liquid He 4 , p n replaces p in equation (A10). The expression is valid provided 
the mean free paths of the rotons and phonons which are assumed to behave like particles 
do not become large enough to violate the assumption of shearing between & layers of 
the fluid (see appendix). Consideration of the phonon and roton mean free paths in liquid 
He 4 indicates equations (All) and (A10) become less applicable to this system as the tem- 
perature is decreased below 1° K. This occurs because (1) near 1° K the phonon mean 
free path becomes larger than the characteristic boundary layer in the liquid and (2) with 
a further decrease in temperature the phonon mean free path becomes large with respect 

to the spacing between the crystal surface and electrode boundary wall. This problem 

3 

does not arise for liquid He in the temperature range of this experiment, 
o 4 

Below 0. 6 K in liquid He , the excitations are nearly all phonons. Furthermore, 
the phonon mean free path becomes considerably larger than the spacing between crystal 
surface and electrode boundary wall. Hence the phonons may be treated as particles in a 
classical Knudsen gas. Here the particles are considered to move from boundary wall to 
the crystal surface and vice versa without colliding with other particles. If it is assumed 
that (1) reflections from all surfaces are entirely diffuse, (2) there is no phonon-phonon 
scattering, and (3) any effect of roton excitation is ignored, it can be shown that the 
crystal decrement is proportional to the phonon density. If the crystal were immersed 
in a classical Knudsen gas, the decrement due to the gas is determined by 
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where 

3 

N number of particles per cm 
S total surface area of crystal 

P classical momentum per particle 

M x mass of crystal 
Fj-^ resonant frequency 

If the classical momentum per particle is replaced by 3p^j i c/4N (ref. 6), the decrement 
due to the phonon gas is 


H = »• 158 > (1 > 

where p ^ is in grams per cubic centimeter and c is the velocity of ordinary sound in 
liquid He^ below 0. 6° K and has the value 238 meters per second. 


APPARATUS 

In the aforementioned setup the torsional crystal and liquid helium samples were 
contained in a vacuum -tight hollow chamber which will be called the crystal chamber 
(see figs. 2 and 3). The walls of this chamber were constructed of copper and were 
surrounded by a thin concentric cylinder of stainless steel. This space between the 
double wall will be referred to as the heat exchange chamber. The heat exchange cham- 
ber was evacuated for some measurements or filled with 10 to 100 microns of helium 

3 

gas pressure for other measurements. The He reservoir was located directly above 

the crystal chamber and was separated from the crystal chamber by a copper partition. 

3 

This reservoir, partially filled with liquid He , was pumped on by a closed pumping 

3 

system consisting of a mercury booster pump, sealed fore pump, and He storage tanks. 
When the primary reservoir was cooled to 0. 95° K, the He^ reservoir could be cooled to 
as low as 0. 38° K depending upon the contents of the crystal chamber. A vapor pressure 
measurement of the He in the reservoir was used to determine the absolute temperature 
of the crystal chamber which was in thermal contact with it. This entire system was 

4 

inserted into a double Dewar system containing liquid He shielded by liquid nitrogen 
(Ng)- The He ^ Dewar functioned as the primary cold reservoir for the He^ reservoir 
and the liquid sample in the crystal chamber. The inner Dewar (shown in fig. 2) 
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Liquid He 4 Dewar- 


Capillary fill tube 
for crystal chamber — 


Primary reservoir 
bath heater- 



Primary He 4 reservoir 
He 3 reservoir pump tube 

Heat exchange chamber 
pump out tube 

Lower radiation trap 
■He 3 vapor pressure tube 


-Vacuum space or heat 
7 exchange chamber 


-He- 5 reservoir pump tube 
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-Copper wall of He 


/ reservoir 


-Copper bottom of He^ 
reservoir and top of 
crystal chamber 


-Crystal in suspension 


-Carbon thermometers 
embeded in crystal 
chamber bottom 

CD-9395 


Figure 2. - Low temperature reservoir system. (Liquid N 2 Dewar has been omitted for clarity. ) 
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containing liquid He was connected to a large vacuum pumping line and other smaller 

lines. This permitted control of the temperature of the primary cold reservoir. The 

3 

heat leak paths available from the primary cold reservoir to the He reservoir can be 

seen in figures 2 and 3. The thin wall stainless steel tubes passing through the top of the 

3 3 

He reservoir represent heat leaks. The larger tube is the support for the He reservoir 

and crystal chamber. A variable heat leak path was provided by regulating the exchange 

gas pressure in the heat exchange chamber from a few millimeters to a good vacuum. 

All electrical wires from the crystal chamber and the heat exchange chamber are carried 

through these tubes to room temperature glass -metal seals. 

The quartz crystal 0. 53 centimeter in diameter and 17. 8 centimeter in length was 

suspended in a crystal holder (fig. 4) . The silver body of the crystal holder served as a 

support for four silver quadrants surrounding but not contacting the crystal. Polytetra- 

fluoroethylene electrically insulates two of the silver quadrants from the silver body. 

The space between the quadrants and the crystal surface (0. 53 mm) is filled with the fluid 

4 3 

of interest, in this case, either liquid He or He . Adjacent quadrants having opposite 
electrical polarity produce electric field gradients in a region about the slots between 
adjacent quadrants. These gradients are responsible for driving the crystal in a tor- 
sional mode. Thus, a capacitative coupling existed between the electrode quadrants and 
the crystal surface permitting alternating current flow. 

The suspension system consisted of two nylon threads contacting the crystal at its 
central strain node on opposite ends of a predetermined diameter. These threads are 
held taut by a suspension spring mechanism. At the tangent points for the two threads, 
two small grooves were sawed into the crystal’s highly polished surface to a depth of 
0. 1 to 0. 2 millimeter. These small grooves provided physical alinement of the crystal 
within the quadrants and also served to aline the electric axes of the crystal lattice with 
respect to the electric field gradients. The required alinement can be attained by adjust- 
ing the slightly offcenter eyelet holes shown in figure 4. The required grooves were 
sawed into the crystal surface by a wire charged with diamond compound, A microscopic 
comparator insured proper positioning and depth of the groove. The groove alinement 
with respect to the lattice axes was achieved by determining the position of the optic axis 
with the aid of a polarized light source and a microscopic slit system. 

Silver caps (see fig. 4) were located about a millimeter from each of the two end 
faces in order to reduce energy dissipation from longitudinal standing waves that may be 
set up with the crystal. This distance is less than one -quarter wavelength of 11 -kilohertz 
sound in liquid helium. 

Utilization of the aforementioned suspension system in a temperature region below 
5° K with the torsional crystal immersed in a few microns of gas provided an oscillator 
having a "Q" of 10 million or higher. Such large Q values offer the possibility of using 
the torsional crystal to establish a highly stable oscillator. 
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Crystal holder body 




Consideration of figure 5 suggests that for a specified crystal voltage the limiting 
sensitivity for the measurement of a constant R x value depends upon the preamplifier 
sensitivity for detection of the bridge output signal and upon the oscillator stability. This 
preamplifier had a gain of 80 decibels and a noise level of less than 1/2 microvolt with 
the input terminal shorted. 

The attainment of a stable oscillator for the driving electromotive force was some- 
what difficult. A very stable oscillator is required for high Q measurements. Experi- 
ments indicate that for moderate values of Q the resonant frequency must be maintained 
within 0. 005 hertz for at least a few seconds so that the bridge can be balanced to yield 
a determination of R within 1/2 percent. This is possible only under ideal conditions 
with a copper shielded beat frequency oscillator having an adequate fine frequency ad- 
justment. Some lower Q data were taken using this instrument. However, high Q data 
required better stability. This stability was obtained by constructing an ultrastable 
oscillator. 

An ultrastable oscillator was constructed by replacing the crystal element in a com- 
mercial frequency standard with a different crystal which was placed in a low temperature 
Dewar system. Because frequency adjustments required inconveniently long times, it 
was later replaced by a commercial frequency synthesizer. The output frequency of the 
synthesizer was selected to be one thousand times greater than the resonant frequency 
of the torsional crystal. This allowed fine frequency control and accurate frequency 
measurements. A divider network was used to reduce the synthesizer frequency by a 
factor of one thousand to obtain the desired crystal resonant frequency. 

In the preceding discussion of bridge output signal stability, it was assumed the tor- 
sional crystal resistance remains constant at resonance while immersed in a liquid at 

4 

constant temperature. This is certainly true for liquid He below the A point. How- 
3 4 

ever, in liquid He or in liquid He above the A point a microscopic bubble, for 
example, will disturb the R value as evidenced by small disturbances observed on the 
null detector. In general less than 1/10 microwatt of power was dissipated by the crystal 
in the sample fluid so that no serious thermal disturbances were observed. 

Temperature of the sample liquid was measured by eight nominal 10-ohm Allen 

Bradley 1/10-watt carbon composition resistors. These resistors were in thermal 

contact with the copper bottom of the crystal chamber. The eight resistors were wired 

electrically in parallel producing a "resistor bank". A maximum of 25 microamperes 

current was passed through this resistor bank. Its electrical resistance was determined 

from potentiometric measurements. Temperature differences of a few tenths of a milli- 

degree could be measured using this resistor bank. The absolute temperature calibra- 

2 

tion for this secondary thermometer was determined from measurements of He vapor 

3 4 

pressure in the He reservoir and under certain conditions from He vapor pressure 

measurements in the heat exchange chamber using octoil-s and mercury manometers. 
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Later the manometers were replaced by fused quartz precision bourdon tube pressure 
gages. These were found to be very satisfactory for vapor pressure determinations 
because of their rapid read out and excellent sensitivity. A graphical fit was then applied 
to these calibration points which determined "resistor bank" resistance as a function of 
temperature. At lower temperatures, a vented germanium thermometer immersed in 
the sample liquid proved very reliable. Temperature cycles from helium temperature 
to room temperature have much less effect on the calibration of germanium thermometers 
than carbon thermometers. 


EXPERIMENTAL PROCEDURE 

The crystal resistance was measured by capacitance coupling to the electrode quad- 
rants thereby eliminating the need for a mechanical contact. The electrode quadrants 
were connected to one arm of a modified Schering capacitance bridge by small electrical 
wires running within the capillary fill tube from the crystal chamber. A capacitance C 

px 

(fig. 1, p. 7) was inserted in parallel with the crystal to keep the arm containing the 

crystal capacitative at frequencies slightly above the resonance frequency. The bridge 

is always used with impedance Z & set equal to impedance Z fo . Hence, when a null 

signal is observed on the oscilloscope, the bridge is balanced and Z = Z . The follow- 

u s 

ing procedure was used to establish the resonant frequency F R and the crystal resistance 
at resonance. 

The bridge was first balanced at frequencies far from resonance where the crystal 
branch containing R , L , and C essentially exhibits an infinite impedance. This was 

X A A 

done by setting R = °° and varying the resultant pure capacitance C in arm S. The 
s s 

capacitance C c at which the bridge balances is then equal to the capacitance C" of the 
crystal arm. Keeping C g = C" results in the bridge being balanced at f r- Under this 
condition a balance at can only be obtained with R g = R x - Hence such a balance 
yields both F R and R x simultaneously. 

The resonant frequency at bridge balance was obtained by viewing a counter and 
timer which has a maximum input frequency rate of 15 megahertz and a maximum gate 
time of 10 seconds. 

The crystal constant K Q was determined by measuring the logarithmic decrement 
directly from the half width of the resonant peak. A plot of R /Z consists ideally of a 
symmetric sharp peak. The half power point of this peak occurs when the bridge resist- 
ance R„ equals 2R . The Q of the peak is defined as 

S X 



6F 2R A 
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where SFg^ is the width of the resonance curve at the half power point. From the 
previous equation and equation (A2) the following equation may be Written: 

R 

K = 5 

7rM6F 2 R 

which is the expression employed in the determination of K Q for the particular temper- 
ature of interest. The previous expression indicates that the frequency must be measured 
with good precision. 

The crystal chamber was evacuated and gradually cooled from room temperature to 

liquid nitrogen temperature over a period of about 36 hours. The primary reservoir was 

o 3 3 

then cooled to about 1. 2 K. Next He gas was condensed into the He reservoir and the 

sample fluid was condensed into the crystal chamber. For liquid He 4 measurements 

considerable care was taken to find a favorable liquid level for reasons considered below. 

4 o 

When the crystal chamber was slowly filled with liquid He at 1. 2 K, the crystal resist- 
ance varied a number of times between maximum and minimum values. A favorable level 
was considered to correspond to a minimum crystal resistance value. This variation in 
crystal resistance at constant temperature with liquid level was observed only in liquid 
He 4 . It was a resonance condition that may be due to second sound waves inducing first 
sound waves in the He 4 vapor at the liquid vapor interface. At certain liquid levels this 
could produce first sound standing waves in the vapor between the liquid surface and the 
crystal chamber top thereby producing additional damping of the crystal. However, if 
the liquid was at a favorable level this unwanted additional loading of the crystal was 
reduced considerably. Not to be confused with the aforementioned problem, unwanted 
first sound loading of the crystal may also be produced by an improper choice of crystal 
circumference for a given crystal length as mentioned previously in the section 
ANALYTICAL CONSIDERATIONS. 


Procedure for Measurements Near the X Point 

4 

After condensation of liquid He was complete in the crystal chamber, the carbon 

3 

resistor thermometer bank was calibrated from He vapor pressure measurements and 

4 

later from He vapor pressure measurements. A series of calibration points was taken 

at nearly constant temperature in the range 1. 9° to 2. 15° K and then later above the X 

point (T^ = 2. 172° K) by correlating resistor bank voltages with the corresponding vapor 

pressure. Over this temperature range below T^ , the liquid He 4 bath reservoirs are at 

3 4 A 

uniform temperature. Both the He and He vapor pressures are rapidly changing 
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functions of temperature so that very accurate absolute temperature determinations are 
possible. Between calibrations, crystal resistance measurements were taken simulta- 
neously with resistance bank voltage measurements. Near 2. 15° K the temperature cali- 
bration of the resistor bank below was stopped. Only measurements of crystal 
resistance against resistor bank voltage were taken as the temperature slowly increased 
through T^. For the small temperature range above the A point, a reliable calibration 
of the resistor bank was obtained by extrapolating from calibration points taken near but 
below the X point. Calibration points taken above the A point were considered less 
reliable because of the poor heat transfer properties of liquid He above T^. The rate of 
change of temperature of the crystal chamber was small enough so that about fifteen 
crystal resistance measurements (and hence V n P n determinations) could be synchronized 
with resistor bank voltages between 2. 15° K and the X point. No data were obtained 
between - 0. 0015° and + 0. 0005° K. The reason for this is that in this small 
temperature interval the ac bridge could not be balanced because the values changed too 
rapidly. Above the X point, measurements were continued to about 2. 20° K. Also here 
it was necessary to depend upon the silver crystal holder body, silver electrode quad- 
rants, and copper wall to produce a uniform temperature in the sample fluid. This ex- 
periment near the X point was performed also with the temperature slowly decreasing 
and with different fluids in the He reservoir and heat exchange chamber with nearly 
identical results. 


Procedure for Measurements in Liquid He^ Below 1° K 

In order to reduce the crystal chamber temperature below 1° K, the He^ reservoir 

was pumped with either a fore pump or a mercury booster pump. The thermometer 

3 

resistor bank was calibrated from He reservoir vapor pressure measurements using a 

McLeod gage and making any applicable thermomolecular correction (ref. 7). Consider- 

4 3 n 

able care was taken to avoid contaminating the pure liquid He with He . At 0. 4 K, 

3 4 

10 parts per million of He in liquid He alters the density of the normal component of 
4 3 

He significantly, since He effectively becomes part of the "normal component. " The 
density of the normal component of liquid He^ at 0. 4° K is only about 0. 5x10“® gram per 

3 

cubic centimeter. Various pump down temperatures for the He reservoir were selected 
by adjusting a metering valve. Measurements of crystal resistance and resistor bank 
voltage were synchronized as the temperature of the crystal chamber was slowly cycled 
from 1. 0° to 0. 38° K. At low temperatures in He^, an accurate determination of crystal 
resistance in "vacuum" is important since it becomes a significant portion of the total 
crystal resistance in liquid He^ (for heat exchange purposes a vacuum was approximated 
by a few microns of exchange gas) . This was accomplished by measuring the crystal 
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4 4 

resistance as a function of temperature in a few microns of He gas before the liquid He 

3 

experiment and in a few microns of He gas after the experiment. The "vacuum" 
resistance for this crystal system varied from 81 ohms at 2. 0° K to 114 ohms at 0. 4° K. 

Procedure for Measurements in Liquid He^ 

3 

The procedure for measurements in liquid He was similar to that already described 

4 3 

for liquid He . As mentioned previously, the liquid level of He in the crystal chamber 

had no detectable effect provided the crystal was completely covered. The crystal 

4 3 

resistance in liquid He at low temperatures was small compared to that for liquid He . 

3 

Therefore the correction for the "vacuum" resistance is less important in He . Liquid 
He is a very poor thermal conductor so that the silver crystal holder body, silver 
electrode quadrants, and copper wall again had to be depended on to provide uniform tem- 
perature for the sample bath. As a consequence, measurements were taken with the 

4 

temperature changing at much slower rates than for liquid He below the A point. 


RESULTS AND DISCUSSION 

4 

Wide divergence of viscosity values for liquid He as a function of temperature was 
reported by Reynolds, et al. (ref. 8). Their figure is shown in figure 6. The wide 
discrepancies are in part due to the very small viscous dissipation that takes place in 
the liquid even though the viscosity coefficient of the normal component is substantial. 
Measurements of such small viscous dissipations are easily disturbed by perturbations 
which would have no effect in a more dissipative liquid. 

The temperature dependence of V n P n near T^ has been measured with better 
precision in this experiment than in previous investigations. The temperature dependence 
of viscosity alone was determined using normal density values calculated by Reynolds, 
et al. (ref. 9). The resultant viscosity values near T^ are in good agreement with 
values quoted in references 10 and 11. These values were obtained with a rotating cylin- 
der viscometer which measures viscosity independent of normal density. The V n P n 
data shown in figure 7 represent two separate experiments covering approximately the 
same temperature range. In one the temperature of the liquid increased slowly, while 
in the other the temperature decreased slowly. These two experiments differed in their 
rates of change of temperature and also in the fluids that occupied the heat exchange 
chamber and the He reservoir. Temperature differences between successive data 
points were measured to better than 0.2 millidegree below T^. The absolute temper- 
ature calibration may be in error by a greater amount. The rate of change of crystal 
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resistance during the measurements together with other considerations indicate it may be 
necessary to shift the data points as much as 1/2 millidegree toward lower temperatures 
on the 58 scale (ref. 12) shown in figure 7 and table I. The viscosity values shown in 
figure 8 and table I were obtained from smoothed V n P n data with the previously men- 
tioned normal density values. Lower values result if the same normal density values are 
used to obtain viscosity values from the V n P n data of both Welber (ref. 13) and Tough, 
McCormick, and Dash (ref. 14). (The normal density values used here are not the same 
as those used to determine viscosity values given in the latter publication. ) Our differ- 
ences in T) n p n values were measured with a precision of better than 0. 15 percent. As 
mentioned previously, data could not be obtained between T^ - 0. 0015° and 


TABLE I. - SMOOTHED VALUES OF 7? n P n AND VISCOSITY ABOUT T x 


Temper- 

ature, 

T 58’ 

°K 

V 

(f*P) 
/or 10" 4 

/'ll’ 3 

(g)/cm 3 
(N) (sec) 

Ag\\ 

\m // 

Density 
of normal 
component, 

Pn’ 

g/cm 3 

(a) 

Viscosity 

V 

pp 

|or ID' 7 (NMsecVm 2 ^ 

IvJ 

2. 195 

3 

82. 

0. 1461 

26. 1 5 

2. 190 

3.78 0 

. 1462 

25. 8 5 

2.185 

3.73 g 

. 1462 

25. 5 ? 

2.180 

3.69 ? 

. 1462 

25. 2 9 

2. 175 

3. 65g 

. 1462 

24. 9 g 

2. 173 

3.63- 

. 1462 

24. 8 g 

2. 172 





2. 1705 

3. 55, 

. 1427 

24. 8 g 

2. 170 

3 

- 53 7 

. 1421 

24. 8 g 

2. 169 

3.47 7 

. 1412 

24. 6 2 

2. 168 

3.42 q 

. 1405 

24. 3 4 

2. 165 

3. 29- 

.1387 

23. 7 6 

2. 160 

3 

•13 7 

. 1360 

23. 0 ? 

2. 155 

3.00 7 

. 1334 

22. 5 4 

2. 150 

2 

,88 6 

. 1310 

22. 0 3 

2.145 

2 

• 77q 

. 1287 

21. 5 2 

2.140 

2 

• 66 1 

. 1265 

21. 0 4 

2. 130 

2 

,46 2 

. 1222 

20. 1 5 

2. 120 

2 

.29. 

. 1180 

19. 4 4 

2.110 

2 

• 15 o 

. 1140 

18. 8 6 

2.100 

2 

.03- 

. 1103 

18. 4 0 

2.090 

1.92- 

. 1069 

18. 0j 

2.080 

1 

o 

CO 

00 

. 1037 

17. 6 5 


a Ref. 9. 
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Figure 8 - Viscosity measurements of He 4 near X point. 


2. 18 


2.20 


23 


+ 0. 0005° K. This 2-millidegree temperature interval where no data were obtained 
seemed to be independent of the rate of change of temperature. A ’’french curve extrapo- 
lation” of the curve passing through the last ten V n P n data points below meets 
at a higher point than a similar extrapolation of values above T^. This discontinuity 
could be due to a very small amount of longitudinal loading exhibited by the torsional crys- 
tal, rather than a viscous property of the fluid. 

The viscosity coefficients determined from V n P n data in this investigation over a 
temperature range extending upward from 0. 95° K are shown in figure 9 and table IT. 

The p n values over this temperature range were obtained from reference 9. The results 

TABLE H. - SMOOTHED VALUES OF 
VISCOSITY FOR LIQUID HELIUM 4 


Temper- 

ature, 

T 58’ 

°K 

Density of 
normal 
component, 

P n’ 3 

g/ciri 

(a) 

Viscosity, 

V 

/ " p \ 
for 10 -7 (N)(sec)/m 2 J 

2.050 

0.0951 

16.6 

2.000 

.0825 

15.1 

1.950 

.0713 

14.0 

1.900 

.0617 

13.4 

1. 850 

.0532 

13.2 

1.800 

.0455 

13.2 

1.750 

.0389 

13.2 

1.700 

.0333 

13.2 

1.650 

.0281 

13.3 

1.600 

.0235 

13. 5 

1.550 

.0194 

13.9 

1.500 

.0160 

14.3 

1.450 

.0133 

14.9 

1.400 

.0108 

15.5 

1.350 

.00883 

16. 2 

1.300 

. 00690 

— 

1.200 

.00419 

— 

1. 100 

. 00220 

23. 0 

1.050 

.00154 

26.9 

1.025 

.00127 

29.4 

1.000 

.00104 

32. 1 

.975 

. 000835 

36.0 

.950 

. 000665 

40.3 


a Ref. 9. 
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Woods and Hollis Hallett (ref. 11) 


? 

5 


Heikkila and Hollis Hallett (ref. 10)J 


^Independent of p n 


Present investigation 

Tough, McCormick, and Dash (ref. 14) 

Welber (ref. 13) 


Smoothed values; p_ values 
from ref. 9 used to deter- 
mine r) n from n n p n values 



reported herein are shown with the values of four other investigations (refs. 10, 11, 13, 
and 14). Except for temperatures near the X point and below 1. 4° K the viscosity values 
are in fair agreement. We do not consider our viscosity values reliable between 1. 1° 
and 1. 35° K because of the previously mentioned problem of the coupling between second 
and first sound at the liquid- vapor interface. The dashed portion of the curve in figure 9 
was obtained by smoothly connecting the data on both sides of this region. Some of the 
measured values were 25 percent higher than the dashed curve. For reasons stated in 
the section ANALYTICAL CONSIDERATIONS, equation (A10) becomes less applicable for 
liquid He^ as the temperature decreases below 1. 0° K. Upon evaluation of the constants 
in equation (A10) appropriate to this experiment, the result is 

A j = 2.333X10" 3 77 n 
1 'n H n 

2 4 

where V n P n is in gram per centimeter per second. The rapidity with which equa- 
tion (A10) becomes less applicable as the temperature is decreased below 1° K is shown 
in figure 10. The values obtained from the equation are compared with the data of Woods 
and Hollis Hallett (ref. 11) which we consider the most reliable data to date. The devi- 
ation between their data and the data in this report which begins near 1° K probably 
occurs because the phonon mean free path becomes larger than the characteristic boundary 
layer depth. In this figure, the normal density values of Bendt, Cowan, and Yarnell 
(ref. 15) are used. The observed peak in this plot at ~0. 6° K was interpreted as the 
point at which the phonon mean free paths become approximately equal to the spacing 
between crystal surface and electrode boundary wall which was 0. 53 millimeter. The 
calculated value of the phonon mean free path at this temperature from Landau's expres- 
sion is about 1.6 millimeters (ref. 1, p. 108). Smoothed values of 

A^xlO 3 
2. 333p n 


are shown in table III (p. 29). 

At temperatures below 0.6° K in this liquid, equation (1) becomes applicable and 
interest is focused on the crystal decrement rather than the square of the decrement. 
Smoothed experimental values of A^ below 0. 6° K are shown in figure 11 and listed in 
table m. In this temperature region there is no longer viscous dissipation in the liquid 
but rather a direct momentum transfer from crystal surface to boundary wall. Here 
momentum is transferred by phonons whose density is highly temperature dependent. 
Equation (1) shows the crystal decrement is proportional to the phonon density; thus a 
determination of the slope of the "line” in figure 11 yields a temperature dependence for 
the phonon density. This temperature dependence is T^’ *, which is in excellent 
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Figure 10. - Measurements on liquid He^. Maximum 
amplitude, ~l/2 micron; maximum velocity, <0.3 
centimeter per second; boundary wall ~l/2 milli- 
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phonon mean free path l p h~0. 07x10“^ millimeter; 
at 0. 95° K: C. B. L. ~0. 013 millimeter, phonon 
mean free path L h ~0.03 millimeter. 
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TABLE m. - SMOOTHED VALUES FOR LIQUID HELIUM 4 


Temper- 

ature, 

T 62’ 

°lc 

(a) 

Phonon density, 
p ph’ 

g/cm 3 

(b) 

A^exp 

Aj exp/Aj cal 

■ 

Temper- 

ature, 

T 62’ 

°K 

(a) 

Density of 
normal 
component, 
P n’ 3 

g/cm J 

(c) 

A^ exp 

A 3 expKlO 3 

2. 333p n 
g/ (cm) (sec) 

0. 375 

0. 3519X10' 6 

2. Ixl0 -7 

0.51 

0.725 

0.044X10 -3 

39. 6X10' 7 


.400 

. 4555 

2.7 

.51 

. 750 

.0630 

44.8 


.425 

.5805 

3.6 

.54 

.775 

.0901 

49.9 

liaxio -6 

.450 

.7296 

4.6 

.54 

.800 

.127 

54.6 

101 

.475 

.9058 

5.7 

. 54 

. 825 

. 173 

58.8 

85.7 

.500 

1.112 

6.9 

.53 

.850 

.229 

62.8 

73.7 

.525 

1.352 

8.5 

. 54 

. 875 

.301 

66.8 

63.4 

.550 

1.628 

10.3 

.54 

.900 

. 391 

70.6 

54.9 

.575 

1.945 

12. 8 


.925 

.500 

74.8 

48.0 

.600 

2.306 

15.9 

— 

. 950 

.632 

79. 1 

d 42.4 

.625 


19.9 


.975 

.790 

83.7 

d 38.0 

.650 


24.7 

— 

1.000 

.974 

88.3 

d 34. 3 

.675 


29.4 

— 

1.025 

1.19 

93.3 

d 3l. 4 

.700 


34.4 

— 

1.050 

1.44 

98.4 

d 28.8 


a Ref. 20. 

b Ref. 1 (p. 68), c = 238 m/sec. 
c Ref. 15. 

d Viscosity coefficient for normal component. (These values will differ from values in table H 
because a different p n set was used). 

agreement with the Landau theoretical expression for the phonon density (ref. 1, p. 68). 
In table HI A^ (experimental) is compared with A^ (calculated from eq. (1)) below 
0. 55° K. The calculated values for A^ might be expected to be higher than the experi- 
mental values since for equation (1) only diffuse behavior was considered. In 
liquid He as the temperature is further decreased below 0. 38 K the crystal decrement 
due to the liquid will become a small fraction of the crystal vacuum decrement. Under 

these circumstances a further decrease in temperature to 0. 1° K will produce no mean- 

3 4 

ingful results. However, if a small amount of He is added to the liquid He , the situ- 
ation will change considerably and the crystal decrement due to the dilute Fermi gas 
will be measurable. 

3 

The results of several "pure" liquid He runs for this investigation are shown in 
figure 12 and table IV. The figure also shows the results of other investigators. The 
continuous curve and dashed curve represent the smoothed values of the corresponding 
investigators whereas the raw data points are being shown from this investigation. 

Below 0. 6° K these data are in excellent agreement with Betts, et al. (ref. 16). Between 
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Viscosity, r?, nP(or 10 _/ (N)(sec)/rrr) 





TABLE IV. - SMOOTHED VALUES OF 


VISCOSITY FOR LIQUID HELIUM 3 


Temper- 

ature, 

T 62’ 

°K 

(a) 

Density, 

P. 

g/cm 3 

(b) 

Viscosity, 

V, 

, pP \ 

/or 10-7 (NHsecjj 

0.360 

0.08214 

65. 1 

.400 

.08216 

57:i 

.450 

.08217 

51.1 

. 500 

.08218 

47.4 

.550 

.08217 

44.6 

.600 

.08216 

42.2 

.700 

.08210 

38.7 

. 800 

. 08202 

36.5 

.900 

.08189 

34.6 

1.000 

.08174 

33.0 

1. 100 

.08156 

31. 3 

1.200 

.08134 

29.2 

1. 300 

.08108 

27.5 

1.400 

.08077 

26.4 

1.600 

.07999 

24.9 

1.800 

.07897 

23.7 

2.000 

. 07765 

22.7 

2.200 

. 07603 

21.9 

2.400 

. 07406 

21.2 

2.600 

.07162 

20.6 


a Ref. 20. 
b Refs. 21 and 22. 


TABLE V. - MEASURED VALUES OF OTHER CRYOGENIC FLUIDS 


Fluid 


VP, 


Density, 

Viscosity, 

Temperature, 


(|iP)(g)/cm 

3 

P. 

V, 

T, 


10“ 4 

(N)(sec) 

Ag\\ 

g/cm 3 

pP 

(or 10' 7 (NKsecJ/m 2 ) 

°K 


m _ 

Lvy 


Ambient air a 

0. 218g 


0. 1164X10" 2 

188- 0 

298±1 

Gaseous helium 4 


. 0512g 


.0622 

CN5 

CO 

77.6 

Liquid hydrogen 13 

21. 0 2 


7. 716 

272. 4 

14.00 

0 

Liquid hydrogen 

19. 8 2 


7.686 

258- 0 

14.00 

0 

Liquid hydrogen 

9-68 9 


7.086 

136. 7 

20.20 


a Passed through anhydrous CaSO^ drying column. 
^Normal hydrogen. 

c Converted to more than 70 percent parahydrogen. 



0. 8° and 1.2° K, these values could be slightly high possibly because of some weak 
resonance condition which is independent of the liquid level in the sample chamber. 

Above 1. 2° K, these values of viscosity fall off faster with increasing temperature than 
the values of Betts, et al. In fact the slope dT 7 /dT of the curve in this region above 

1. 2 degrees is in agreement with that of Zinov'eva (ref. 17) and Osborne, et al. (ref. 18). 

The values lie between those of Betts, etal. and Zinov'eva. If these values were in 

error, they would be expected to be high. The reason for this is that any significant 

error would probably be the result of a nontorsional component of crystal motion. Such 

motion increases the energy dissipated by the crystal thus adding to the decrement caused 

by pure torsional motion. For this reason we expect these values above 1. 6° K to be 

more acceptable than those of Betts, et al. The experimental precision for rjp of liquid 
3 

He is about 0. 2 percent. 

In addition, the same technique was used to measure the viscosity-density product 
for other fluids. These values are presented in table V to provide an opportunity to 
compare this method with other techniques. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 25, 1967, 

129-01-02-04-22. 
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APPENDIX - DEVELOPMENT OF EQUATIONS OF INTEREST 

The logarithmic decrement of a piezoelectric crystal or of any mechanical vibratory 
system is defined by the quotient of energy dissipated per cycle W d divided by twice the 
energy of vibration W . In this case the energy of vibration is simply the maximum 
kinetic energy attained by the crystal during a cycle. Thus, 


2W v 

where A is used to denote the logarithmic decrement. A fundamental piezoelectric 
equation relates the logarithmic decrement to the electrical resistance and frequency of 
the crystal at resonance (ref. 5): 


R x = K o MF R A 


(A2) 


where R x is the crystal resistance at resonance in ohms, K q at a specific temperature 
is dependent upon the crystal, fluid dielectric properties, and electrode geometry, M is 
the crystal mass in grams, and Fj^ is the crystal resonant frequency. K Q is temper- 
ature dependent and generally decreases with decreasing temperature; however, changes 
in K q were found to be less than 1/2 percent between 4. 2° and 0. 4° K. 

The energy dissipated per cycle is obtained from the power loss across the crystal 

resistance R at resonance 
x 


W d = 



(A3) 


where W d is in joules per hertz if I is the root mean square current in amperes. The 
energy of vibration of this crystal may be calculated from the maximum kinetic energy of 
the torsional rod or crystal, that is, 


W v = \ MV M 


(A4) 


where V 


is the maximum velocity of a point on the rim of the cylinder. Substituting 


results in the following expression for the maximum velocity in centimeters per second 
of a point on the rim of the torsional crystal: 

2EKy 2 10 7//2 

V=V M = - 


R 


33 



where E is the root mean square voltage across the crystal at resonance (E = IR) in 
volts. 

If the maximum velocity of a point on the rim of the cylindrical crystal is known, 
the amplitude of this point or any such arbitrary point on the rim may be determined. 
The instantaneous angular velocity of any point on the rim is given by 


M 

6 = — — cos wt 


where r Q is the radius of the crystal cylinder and, hence, 

y 

6 = — M sin cot 

r co 
o 


Because co = 2 ttFj^, the maximum linear displacement in centimeters for an arbitrary 
rim point will be 

In order to determine rjp in terms of measured quantities, it is convenient to find 
W y in terms of 77 . A relation between the viscosity coefficient and the crystal system 
decrement may be developed after first considering the motion and forces on the fluid 
adjacent to the surface of a body moving with respect to the fluid. It is assumed here 
that a viscous laminar flow is being dealt with which simply means that each 6 layer of 
fluid moves along with the moving boundary surface and exerts a shearing stress or 
tangential force on its neighboring 6 layers. The term 6 layer here refers to a layer 
of fluid having an arbitrarily small thickness normal to the direction of motion of this 
layer. For the case of laminar flow, the tangential force per unit area is proportional 
to the velocity gradient normal to the 6 layer. The proportionality constant is called 
the viscosity coefficient, that is, 


shearing stress 

V = 

velocity gradient normal to flow 

Consideration of mass conservation and the rate of change of linear momentum for 
an element of fluid provides a fundamental set of equations known as the Stokes -Navier 
equations. These equations provide a convenient starting point for nearly all hydro- 
dynamic problems. The problem being solved herein has cylindrical geometry. How- 
ever, the experimental conditions are such that the motion of a fluid caused by an oscil- 
lating plane affords a good approximation to the actual motion. First, consider a 
boundary value problem consisting of an oscillating plane in an infinite fluid. The 
Stokes -Navier equations in rectangular coordinates for an incompressible, constant 
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density fluid with zero body forces are (ref. 19) the following: 

P 


Du _ 3p , Jd 2 u 3 2 u h 3 2 u \ 


Dt 


0X 


Dv = _ 9p + r,/ 


3 v , 3 v . 3 v 


Dt 


3y 


(3x 


„ Dw _ 3p . „(3 2 w , 3 2 w , 3 2 w\ 
p + VI + — - + 


Dt 


3z 


\3x 



(A 5a) 


(A 5b) 


(A5c) 


where u, v, and w are the x, y, and z components of the fluid velocity, p is the 
pressure, and 


_D___3_ + u3 + v3 + w3 
Dt 3t 3x 3y 3z 


Next determine the motion and forces of the fluid adjacent to the surface of an infinite 
plane oscillating with respect to the fluid. Consider the oscillating plane as lying in the 
y-z plane of this coordinate system and oscillating in the y-direction. The fluid lies to 
the right of this boundary plane along the x-axis as shown in the following sketch: 


y 



Now apply the Stokes -Navier equations to this problem. As set up, within the fluid 
u = w = 0 and conditions at a point depend only on the x-coordinate of that point. Equa- 
tion (A5c) may be eliminated since the quantities are all zero. Equation (A5a) reduces to 
the simple relation 3p/3x = 0, while equation (A5b) becomes 


P 


— = t] d v 


3t 


dx“ 


The boundary condition at x = 0 is based on the assumption that the first layer of fluid 
atoms immediately adjacent to the surface of the plane moves with the plane and exhibits 
no slippage, that is, at x = 0, v = V m cos cot where V m is the maximum velocity of the 
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plane boundary. Also due to friction, v = 0 at x = 

The aforementioned differential equation is identical in form to certain heat flow 
and vibrating string problems and may be handled by separating variables. The solution 
may be written as 

v(x, t) = V m e“''' x cos (cot - Ax) (A6) 

where v(x, t) specifies the velocity of a particle in the fluid at a distance x from the 
plane surface at time t and 



Hence, the velocity of fluid particles is reduced exponentially as the normal distance 
from the infinite plane increases. At a distance x = 1/A the velocity is reduced to one 
e^ of its value at x = 0. This distance, 1/A, is the thickness of the Characteristic 
Boundary Layer and henceforth will be referred to as C. B. L. 

An expression for the velocity of a given fluid particle as a function of its distance 
x from the plane surface has been obtained. Under certain conditions, a cylindrical 
surface may be treated hydrodynamically as though it were comprised of a large number 
of small plane areas. These conditions are the following: 

(1) The C. B. L. is very small in comparison with the radius of the cylinder r . 

(2) The maximum amplitude of any point on the surface of the cylinder is very small 
in comparison with r . 

(3) The length of the cylinder 2L is large in comparison with the radius r . 

This approximation may be further extended to the case of a cylindrical rod in tor- 
sion provided that the following conditions are met: 

(4) The nonzero shearing stress due to 3v/3z is very small compared with the 
stress due to the velocity gradient 3v/3n normal to the cylindrical surface. 

(5) The nonzero shearing stress on the end faces of the torsional rod or crystal due 
to 3v/3R 0 is very small as compared to the stress due to the velocity gradient 3v/3n 
normal to the end faces, where R g is the radial coordinate of the end faces. Note that, 
for the crystal being considered, v = 0 at R = 0 and v = V M cos cut at R = r ; 
thus, 3v/3R g / 0. 

In general, conditions (1) to (5) are not assured; but for the case of a piezoelectric 
torsional crystal having a resonant frequency of 11 kilohertz, these conditions are 
certainly satisfied for the most part. One would suspect this to be the case since high 
frequency implies a very small C. B. L. which in turn implies enormous velocity 
gradients normal to the crystal surface as contrasted to the other gradients 3v/3z and 
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3v/3R g . The C. B. L. of this crystal in liquid He^ was only 2. 6xl0 - ^ centimeter at 
1. 30° K. On the surface of the cylinder, the maximum values of the various velocity 
gradients are 0v/3n) max = Vj/, 0v/3z) max = (tt/ 2L)V m , and 0v/3R e ) max = 8^ 

The velocity gradient normal to each element of area of this crystal surface may be 
obtained by differentiating equation (A6) with respect to n (which has replaced x) and 
noting that the velocity V m on the surface of the cylinder is a function of position. 

Essentially the variable x of the oscillating plane problem has been replaced by the 
variable n where n is measured everywhere normal to the cylindrical surface including 
the ends and has a value n = 0 at the crystal surface. Consider a torsional rod of length 
2L which is oriented such that the longitudinal axis of the rod lies on the z -coordinate 
axis with the end faces of the rod at z = L and z = -L (as defined before, V M is the 
maximum velocity of a point on the rim). The torsional motion for an undamped rod is 
described by the following equation: 


d 2 e _ l d 2 e 
3z 2 v 2 at 2 


where 6 is the polar angle and Vg is the velocity of shear in quartz. The solution of 
this equation corresponding to the lowest fundamental mode with strain node at z = 0 is 


9 = 8 


M 




and 



This mode of vibration is the same one which was previously called the fundamental mode 
of vibration possessing the frequency Fp given in equation (All). The angular velocity 
co may be written as 


7TV- 

co = 27rF„ = — £ 
R 2L 


Provided that conditions (1) to (5) are met, the fundamental definition of viscosity enables 
one to formulate an expression for the energy dissipated per cycle for an element of 
area, either on the cylinder or an end face. Thus from fundamental considerations the 
following may be written: 
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dS dt 


r2tr/w 

dW d = 71 I ' 


v(z,R e ,n=0,t) /— 
\ 3n> 


n=0 


Now using equation (A6) and noting that V m is a function of position, that is, 


v = V m (z) e"^ n cos (cut - An) 


for the cylindrical surface or 

v = V m (R e ) e - '*' 11 cos (cut - An) 

for the end face result in the velocity gradient for an element of area on the crystal 
surface 


— ) = V (sin cut - cos cut) 

9n /n=0 


where for the cylindrical surface 


m 


= V z > = V M 



and for the end faces 


V m = V m «e> = Vm 


Then the energy dissipated per cycle for the entire crystal surface including the ends is 
’L pt=2ir /(. u 

W d = -277A ' 


I 


Vm sin (iZ/ cos wt ) 


V M sin ( (sin cut - cos cut) 


,2L) 


2ttt q dz dt 


r T 0 

-2r?A / / R g 0 M (cos cut)R e ^ M (sin cut - cos cut)27rR e dR g dt (A8) 

Jo Jo 


1 /2 

where A = (pcu/2p) ' , 1/A isC. B. L. , is the maximum velocity of a point on the 
rim of the cylinder and is constant for a given driving voltage across the crystal, (? M 


38 



is the maximum angular velocity corresponding to a point on the end face or rim of the 
cylinder and is also a constant, and R g is the radial distance to an arbitrary area strip 
on the end face and varies from 0 to r Q (cylinder radius). Integration of equation (A8) 
yields 


W d = 



(A9) 


where S is the entire surface area of the crystal including the two ends. Substituting 
equation (A7) in equation (A9) and then equations (A9) and (A4) in equation (Al) result in 


A 2 = *S 2 VP 
l p 

F R M 


(A10) 


Upon further substitution from equation (A2), the following is obtained: 


VP = 


< R x - R x, vac> 2 
K o s 2 f r’ 


(All) 
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